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We calculate the light quark spectrum of lattice QCD in the quenched approximation using staggered 
quarks. We take the light quark mass, infinite volume, continuum limit. With non-linear chiral 
extrapolations, we find that the nucleon to p mass ratio is mN /rrip = 1.254 ± 0.018 ± 0.028, where 
the errors are statistical and systematic (within the quenched approximation), respectively. Since 
the experimental value is 1.22, our results indicate that the error due to quenching is ^5%. 
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The calculation of the light hadron spectrum is one of 
the main goals of lattice QCD. However, the computed 
value of the nucleon to p mass ratio has persistently been 
too large. This has been true for calculations done with 
dynamical fermions, or in the quenched approximation 
101, and for both the Wilson and staggered methods of 
putting quarks on the lattice. In such calculations, 
one must use a finite-volume box and a non-zero lattice 
spacing. In addition, because of the iterative algorithms 
used to calculate the quark propagators, quark masses 
larger than those of the up and down quarks are nec- 
essary. Thus, before comparing with experiment, one 
must take the infinite volume, zero lattice spacing, and 
light quark limits to remove systematic errors. Our ex- 
tensive series of calculations with staggered quarks in 
the quenched approximation provides good control over 
each of these extrapolations [^,|| . Here we present a new 
analysis of the chiral extrapolation fully consistent with 
quenched chiral perturbation theory |^ . 

It is important to demonstrate that calculations with 
Wilson and staggered quarks have the same continuum 
limit. Butler et al. ^ calculated the quenched Wilson 
quark spectrum using three lattice spacings. Their re- 
sults had a statistical accuracy of 5-6%, and they claimed 
good agreement with the real world for a number of 
hadron mass ratios. However, there are concerns about 
how well they were able to control their extrapolations 
|0,Q. Recently, the CP-PACS collaboration has pre- 
sented preliminary results with Wilson quarks, using four 
values of the gauge coupling on large lattices [||. 

In Table |, we list the coupling, volume and size of 
each or our ensembles of lattices. (Computational details 
are found in Ref. B.) We used five quark masses for 



each ensemble, each a factor of two different from the 
next heavier or lighter. To increase our statistics, we 
calculated propagators from every eighth time slice, so we 
have four to eight sets of hadron propagators per lattice 
(depending on Nt). These propagators are correlated, so 
they are blocked together for further analysis. 

Now let us turn to the sources of systematic errors. 
The initial motivation for these calculations was to un- 
derstand finite volume effects. Q. We have carried out 
detailed studies at gauge coupling 6/.g^ = 5.7 using lat- 
tices with spatial dimensions Ns — 8, 12, 16, 20, 24, 
with the same five quark masses on all five lattices. At 
6/g2 = 5.85, we use three sizes, Ns = 12, 20 and 24. The 
last column of Table | contains the box size in fm based 
on a scale set by the p mass extrapolated to the physi- 
cal quark mass. In each case, the smallest size is 1.8 fm, 
considerably smaller than the 2.7 fm size for &/ = 6.15. 
We expect that in a small volume the mass will increase 
as the hadron becomes squeezed by the box, and that 
this effect increases for smaller quark masses. For the p 
there is very little evidence for finite size effects. Com- 
paring the smallest and largest sizes for the three lightest 
quark masses, the largest difference is 2.5 ± 1.3% and the 
p is lighter on the smaller size. We expect somewhat 
larger finite size effects for the nucleon. For Qj — 5.85 
we find a 3.2 ± 1.1% effect for the lightest quark mass; 
however, at 5.7 for the corresponding quark mass the 
effect is only 1.0 ± 1.2%. For small box sizes, there is 
evidence that finite-size effects fall like \/V |jl^, where 
V is the spatial lattice volume. For larger boxes, an ex- 
ponential decay proportional to exp(— m,r-^^) is expected 
pl| . By the former consideration, the finite size effect 
for our 2.7 fm = 6.15 lattice should be smaller by 
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a factor of (2.7/1.8)^ « 3.4. By the latter, there should 
be a decrease by a factor of « 5. Thus, we expect finite 
size effects smaller than 1% for 6/.g^ = 6.15, and even 
smaller effects for the two largest volumes at 6/(7^ — 5.7 
and 5.85. 

The chiral extrapolation requires the most care and 
controls the final error. It is based on a small mass ex- 
pansion, but since we cannot generate accurate results 
for small quark masses, we are forced to use intermediate 
and heavy quark masses. In this mass region the contri- 
bution of the higher order terms in the chiral expansion 
can be significant, and the extrapolation to the physical 
light quark masses has to be performed very carefully. 
Another complication arises from the quenched approx- 
imation. Quenched chiral perturbation theory {QxPT) 
1^ differs from the ordinary chiral perturbation theory 
(xPT) that applies to the real world. 

Our lightest three masses cover a range of a factor of 
four, and cannot be fit to a linear function with even 
a marginally acceptable confidence level (CL). (The chi- 
ral fits are done using the full covariance matrix of the 
hadron masses for different quark masses.) We began 
by trying a dozen different fitting forms (Table ||) . The 
term proportional to mq appears at tree-level (in both 
XPT and QxPT) and is therefore expected to be the 
most important correction to the chiral limit for a "rea- 
sonable" range of masses. The other niq dependent terms 

should be thought of as one-loop corrections: m^^^, m^, 
and rriq log niq are standard terms which appear in both 
xPT and QxPT, while the more singular (at ruq — 0) 
terms proportional to , 



q and TTiq log niq arise only in 
QxPT, and are due to 77' loops However, we reject 
fits 1,2,3,4 and 10, which include an ^yruq term, because 



the coefficient of y/mq found by the fits is opposite in 
sign to what is expected from QxPT. In Fig. 1, we show 
seven fits for the nucleon at — 6.15 that have good 
CL. As can be seen in the inset, several fits (2, 3, 4, 10) 
decrease sharply at small quark mass. However, a 
term with the sign implied by QxPT would cause a sharp 
increase. 

In addition, due to flavor breaking, the ^rn^ term is 
not really appropriate to quenched staggered quarks. It 
is actually the flavor singlet pion that appears in 
QxPT, and its mass is not proportional to y/niq at fi- 
nite lattice spacing. Thus, we should use instead a term 
proportional to the mass of the non-Goldstone pion, com- 
monly denoted TOttj • Since QxPT gives us only a rough 
idea of the values of the proportionality constants, we 
consider a range of values. In practice, we define, for 
fixed Ai and A2, 
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where "phys" stands for the physical values and the other 
quantities are values computed at given quark mass and 



lattice spacing. We then fit 
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and 



to functions 8 



and 12 (Table ||) for various values of Ai and A2 obeying 
0.0 < Ai < A2 < 0.4, which is the expected range of 
values from QxPT 

Table HI contains the combined CL of these chiral fits 



for our five lattices with the weakest couplings and largest 
volumes: Ns = 32 at G/g^ = 6.15, and Ns = 20 and 24 
at 6/5^ = 5.7 and 5.85. Although many of these fits 
are marginally acceptable, none of the CL are very good. 
This may be due to the fact that QxPT to order 
would require all the terms that appear in either Fit 8 or 
12; while we are limited to a maximum of 4 fit parameters 
since we have only 5 masses at our disposal. It is never- 
theless encouraging that almost all the fit parameters are 
of the rough size (within a factor of 2) and sign predicted 
by QxPT |l|Jll- The exception is the term in the 
p fits: its coefficient is more than an order of magnitude 
smaller than the size suggested in Such agreement 
is as good as could be expected, since the QxPT pre- 
dictions are currently based either on rather arbitrary 
guesses of parameters {e.g. the n^J"^ case for the p) or 
on parameter estimates taken from real- world (i.e., not 
quenched) data and determined only up to large errors. 
As expected, the linear term in niq is dominant; all other 
terms give small corrections except at the largest values 
of the quark mass. The fit parameters are relatively sta- 
ble: As the fit (8 or 12) or A values are changed, the 
parameters change by at most ±25% from their average 
values, while the masses extrapolated to the chiral limit 



are considerably more stable (see Table IV). The range 
of results for the acceptable fits gives us an estimate of 
the chiral extrapolation error. All the other fit functions 
in Table || are either inconsistent with QxPT or have 
very low CL: the highest of these is ~ 0.0001. 

The final extrapolation is in the lattice spacing a. Al- 
though we expect O(a^) errors in general for the stag- 
gered action, flavor symmetry breaking gives m^^ ~ -^'^^ 
at niq — [ p5[ , which implies that contributions to 
the nucleon and p could produce 0{a) terms in /nip. 
Since, however, we are attempting to remove the to^j 
dependence, we believe it is still reasonable to fit to the 
quadratic form C + B{amp)^ for our central values. In 
Fig. H, we show the lattice spacing extrapolation based 
on chiral fit 8 with Ai — 0.0 for the nucleon and fit 12 
with A2 = 0.1 for the p, which give the highest CL. 
For the tt chiral extrapolation, which is needed to de- 
termine the physical light quark mass, we use the form 



mi = am„ -|- bmi + cniq 



dniq In niq. When we keep 



only the three weakest couplings, the continuum extrap- 
olation with these parameters has a good CL (0.74) and 
gives a result which is close to the overall average of the 
continuum extrapolations of all the chiral fits. For each 
of the two intermediate couplings, we plot the two largest 
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volumes, but only the largest volume was included in the 
fit to the a dependence. Extrapolating to the continuum 
limit we get 1.251 ± 0.018. 

To obtain our central value, we look at all (quadratic) 
continuum extrapolations of the three weakest couplings 
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coming from any of the chiral fits in Table III with CL 
greater than 0.04. (The continuum extrapolations all 
have CL > 0.34.) Figure || shows the results of these 
extrapolations of /m'p as a function of the CL of the 
continuum extrapolation. After averaging over all fits 
in Fig. |, we obtain m^/mp = 1.254 ± 0.018 ± 0.022, 
where the last error is the standard deviation over the 
fits, which we take as a combination of chiral extrapola- 
tion and continuum extrapolation errors. The result is 
rather insensitive to the cut on chiral CL: changing the 
cut to 0.005 or 0.06 (a cut greater than 0.06 would rule 
out all the p fits in Table [II) changes the central value 
by 0.010 or 0.008, respectively. 

To explore further the error in the continuum extrapo- 
lation, we then include the strongest coupling point and 
repeat the analysis. Averaging over all fits with chiral 
and continuum CL greater than 0.04 gives 1.248±0.016± 
0.008. We then change the continuum extrapolation to 
include the higher power (amp)**, as well as (amp)^. By 
averaging as before, we obtain 1.266 ± 0.020 ± 0.021. 

Because the true values of Ai^2 are not known, we can- 
not definitively subtract off the to^^ dependence. There- 
fore, there may remain a small 0(a) term in m'j^/m'p. To 
study this effect, we add a linear term in a to our cen- 
tral continuum extrapolation fits above. (These are now 
constrained fits.) The coefficient of the linear term is al- 
ways consistent with 0, with the errors giving a bound on 
its magnitude of the size expected from a small residual 
m,r2 contribution. Averaging the extrapolated values of 
m'^/m'p gives 1.269 ± 0.096 ± 0.016. 

These considerations lead us to include additional er- 
rors of 0.015 (the effect changing the continuum extrap- 
olation) and 0.010 (the effect of changing the cutoff on 
the CL). Combining them in quadrature with the 0.022 
determined above, gives a total systematic error due to 
chiral and continuum extrapolations of 0.028. 

In summary, the chiral extrapolation is the most del- 
icate issue in our computation of the nucleon to the p 
mass ratio. The simple linear chiral extrapolation is 
ruled out. Our results are reasonably well described by 
fits motivated by QxPT. However, the CL of such fits 
is a slowly varying function of the parameters Ai and 
A2, which therefore are not determined in our procedure. 
Instead, we average over reasonable ranges of these pa- 
rameters. Fortunately, the continuum values of mj^/mp 
produced are not strongly dependent on the parameters. 
Taking into account the variance over the fits, we ob- 
tain TOjv/TOp = 1.254 ± 0.018 ± 0.028. Comparison with 
the observed value of 1.22 indicates that the effects of 
quenching are less than about 5% at the 1 a level. 

This work was supported by the U.S. DOE and NSF. 
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FIG. 1. The nucleoli mass vs. quark mass with = 6.15. 
We plot fits 12, 8, 9, 10, 3, 4 and 2 (listed in order of decreasing 
chiral limit), all of which have CL greater than 0.45. 
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FIG. 2. niN/mp at the physical quark mass vs. {airip)^. 
The points come from fit 8 for the nucleon and 12 for the 
p, with Ai = 0.0, A2 — 0.1. The horizontal line indicates 
the physical value. The solid straight line shows a linear (in 
(arrip)^) extrapolation that includes only the three smallest 
lattice spacings; the dot-dashed straight line, all four lattice 
spacings. The curve is a higher order fit to all four spacings. 
The extrapolated values near arup = have been spread hor- 
izontally for clarity. 



FIG. 3. rn'j^/m'p extrapolated to a = vs. the CL of the 
fit used for this extrapolation. The points represent different 
choices for the chiral fits of the p and nucleon and different 
values of the parameters Ai and A2. All extrapolations in a 
are quadratic and keep only the three weakest couplings. The 
combined CL over these couplings of each chiral fit included 
here is greater than 0.04. 



TABLE I. Lattices used for spectrum calculation 



6/5' 




size 


# aruq 


aN, (fm) 


5.54 




163 


X 32 


205 0.02-0.32 


4.96 


5.7 




8^ 


X 48 


605 0.01-0.16 


1.83 


5.7 




12^ 


X 48 


405 0.01-0.16 


2.75 


5.7 




16^ 


X 48 


405 0.01-0.16 


3.66 


5.7 




20^ 


X 48 


205 0.01-0.16 


4.58 


5.7 




24^ 


X 48 


199 0.01-0.16 


5.49 


5.85 




12^ 


X 48 


205 0.01-0.16 


1.76 


5.85 




20^ 


X 48 


205 0.01-0.16 


2.93 


5.85 




24^ 


X 48 


200 0.01-0.16 


3.52 


6.15 




32^ 


X 64 


115 0.005-0.08 


2.65 








TABLE II. 


Our fitting functions 




Fit 1 








M + arnj'^ 




Fit 2 








M + arriq + hruq 




Fit 3 








M + mn'J'^ + hniq 


3/2 

+ cmq 


Fit 4 








M + arn^^ + hm,q 


+ cml 


Fit 5 








M + amq 




Fit 6 








M + aniq + hrr?q'^ 




Fit 7 








M + aniq + 6mg 




Fit 8 








M + arUq + brriq^^ 


+ cml 


Fit 9 








M + arriq + bruq lo 


gruq 


Fit 10: 






M + amriq^^ + brUq 


+ CrUq log TUq 


Fit 11: 






M + arUq + 6mg lo 


grUq 


Fit 12: 






M + arUq + hml + crriq log niq 



TABLE III. Combined CL of cliiral fits 



Fit 


A = 0.0 


A = 0.1 A = 0.2 


A = 0.3 


A = 0.4 


Nucleoli Jackknife fits 


8 


0.186 


0.140 0.103 


0.075 


0.054 


12 


0.086 


0.050 0.028 


0.015 


0.008 


Rho Jackknife fits 


8 


0.048 


0.041 0.035 


0.031 


0.028 


12 


0.059 


0.060 0.054 


0.046 


0.037 



TABLE IV. Dependence of nucleon and p masses on chi- 
ral fit at &/g^ = 6.15. mjv = niN + Xim^^ extrapolated to 

niq = 0. "niN" = 'rfiN — Aim^j; whore rn,r2 is separately 
extrapolated to = with the same fit. Similarly for p. 





Fit 


Ai = 0.0 


Ai = 0.2 


Ai = 0.4 




8 


0.404(6) 


0.427(5) 


0.450(5) 


TUN 


8 


0.404(6) 


0.407(5) 


0.410(5) 


rfiN 


12 


0.409(5) 


0.433(5) 


0.457(5) 




12 


0.409(5) 


0.413(5) 


0.417(5) 


rhp 


8 


0.320(4) 


0.343(4) 


0.366(4) 




8 


0.320(4) 


0.323(4) 


0.326(4) 


rhp 


12 


0.319(3) 


0.343(3) 


0.367(3) 


"rup" 


12 


0.319(3) 


0.323(3) 


0.327(3) 



5 



